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! Abstract 

o . 

£f~) ■ Using the SANC system we study the one-loop electroweak standard model pre- 

dictions, including virtual and real photon emission, for the decays of the on-shell 
— ' vector boson, W — > ^(7). The complete one- loop corrections and exact photon 

^p" 1 ' emission matrix element are taken into account. For the phase-space integration, 

Qh! the Monte Carlo technique is used. This provides a useful element, first for the 

I evaluation of the theoretical uncertainty of PHOTOS. Later we analyse the source of 

the differences between SANC and PHOTOS and we calculate the additional weight, 
• i-H . which once installed, improves predictions of PHOTOS simulations. We can conclude 

that, after the correction of the weight is implemented, the theoretical uncertainty 
C$ ' of PHOTOS simulations due to an incomplete first-order matrix element is reduced 

to below 2f, for observables not tagging the photon in a direct way, and to 10% 
otherwise. This is interesting for applications in the phenomenology of the ongoing 
LEP2 and future LC and LHC experimental studies. 
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1 Introduction 



The SANC project of Ref. pQ has several purposes. The intermediate goal is to summarize 
and consolidate the effort of the last three decades in calculating Standard Model radiative 
corrections for LEP, in a well organized calculational environment, for future reference. 
However, it is aimed not only at training young researchers and students, but at some 
remaining calculational projects for LEP as well. Recently it was used for that purpose 
in 0. 

An important lesson from the LEP experiments [3 is that the desirable way of provid- 
ing theoretical predictions is in the form of Monte Carlo event generators. This aspect has 
been taken into account in the development of SANC from an early stage of its development. 

The currently available version can construct one-loop spin amplitudes for the decays 
of the gauge bosons W and Z and of the Higgs boson H. For the moment, SANC features 
single real-photon emission, in the calculations of the total rate and decay spectra of the 
B — > //(7) process. The complete spin polarization density matrix of the decaying boson 
is taken into account as well. 

The integration, with the Monte Carlo method, over the three- (two)-body final state 
is done without any approximation (in particular, the small-mass approximation is not 
used). The program provides MC events with constant weight (unweighted events). The 
whole system is, therefore, fairly self-contained and complete. 

It is of the utmost importance for such a system to reproduce known results precisely; 
this step of its development was completed in ref. Let us turn now to its application: 
first find an approximation of the relatively large terms missing in PHOTOS [4] for generation 
of bremsstrahlung in W decay, later install them into this program, and finally verify, again 
with SANC, that the PHOTOS algorithm for the generation of bremmstrahlung in decays, 
after modifications, works indeed better. 

Our paper is organized as follows: in the next section we recall, from the set of 
observables chosen for tests and we explain the input parameters used in SANC and PHOTOS. 
Section 3 is devoted to a discussion of the SANC matrix elements for QED bremsstrahlung. 
The gauge-invariant part of this matrix element, which is expected to be responsible for 
the PHOTOS discrepancies with exact matrix element, is removed from SANC first. These 
terms are approximated by simple formula and installed to a PHOTOS. Section 4 describes 
comparisons of full SANC with PHOTOS with the new correcting weight introduced. A 
summary, section 5, closes the paper. 

2 Initialization set-ups for SANC and PHOTOS runs 

In the following sections we compare predictions from the programs SANC and PHOTOS. It 
is essential that the initialization be identical in all cases and close to the physical reality; 
in particular the following options are set in the two programs: 

• In SANC we switch off the EW part of the radiative corrections. The soft/hard 
photon limit is kept at 0.005 of the decaying particle mass. 
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• In PHOTOS we switch off the double bremsstrahlung corrections. The soft/hard 
photon limit is kept at 0.005 of the decaying particle mass. For the generation of 
the Born-level two-body decays, we use the Monte Carlo generation from SANC. 

To visualize the differences (or the agreement) between the calculations, we choose 
a certain class of (pseudo-)observables, more precisely the one-dimensional distributions, 
which are quite similar to the ones used in the first tests of PHOTOS reported in [SJ. To 
visualize the usually small differences, we plot ratios of the predictions from the two 
programs rather than the distributions themselves. 

List of observables: 

• -A- Photon energy in the decaying particle rest frame: this observable is sensitive 
mainly to the leading-log (i.e. collinear) non-infrared (i.e. not soft) component of 
the distributions. 

• -B- Energy of the final-state charged particle: as the previous one, this observable 
is sensitive mainly to the leading-log (i.e. collinear) non-infrared (i.e. not soft) 
component of the distributions. 

• -C- Angle of the photon with final-state charged particle: this observable is sensitive 
mainly to the non-collinear (i.e. non-leading-log) but soft (i.e. infrared) component 
of the distributions. 

• -D- Acollinearity angle of the final-state charged particles: this observable is sen- 
sitive mainly to the non-collinear (i.e. non-leading-log) and non-soft (i.e. non- 
infrared) component of the distributions. 

3 Complete and truncated matrix elements of SANC 

The Feynman diagrams for the decay 

W ± (Q, A) -> l(p h A/) + u(p u , A,) + 7 (fc, a) (1) 

are shown in fig. 1 (unitary gauge), where (pi, A/) and (p v , \ u ) denote the four-momentum 
and helicity of the fermion I and neutrino, respectively, (k, a) is the four-momentum and 
helicity of the photon, and (Q, A) the four- momentum and helicity of the W boson. 
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Figure 1: Feynman diagrams 



In the notation of ref . jH] , the corresponding parts of spin amplitudes take the form: 



Q, 



BxuxAPhQiPx 



Qi 



2k- Pl 



Pi 



p=± 



Q 



+ t^ jP (p,, m,, k, 0, fc, 0)5^ (A;, Q, p v ) 
where we have introduced the following notation : 



(2) 



B\ 1 ,\ 2 (pi,Q,P2) = ^=u(pi,Ai)e^(Q)(l+7 6 )u(p2,A 

^Ai,A 2 (Pi' m i' fc ' ^,P2,m 2 ) = u(p 1 ,\ 1 )Z f (k)u(p 2 ,\ 2 ) , 
Sx 1 \ 2 b a (k,p) = UZ lM (p,m,k,0,p,m) , 



2j , 



(3) 



and Qiy are the electric charges of the fermion I and the W boson, respectively, in 
units of the positron charge, e^(k) and e^(Q) denote respectively the polarization vectors 
of the photon and the W boson. 

This amplitude can be divided into three gauge-invariant parts: 

• Infrared-divergent part: 



MZ M {k,Q,p v ,pi) 



to 



Qi 

2k- Pl 



b a {k,Pi) 



(4) 



W 



Q 

2k-Q 



(b a (k,pi) + b a (k,p v )) 



B l l>Xl ,iPhQ,Pu) 
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• Collinear-divergent part: 

Ml Xv>Xi (k,Q,p v , Pl )® = -^Y,UlM m i^ k ^)K-K^Q^) ( 5 ) 

/ K ' Pl P =± 

• Finite part : 

+ UZ up (p h m h k, 0, fc, 0)5^ (k, Q, pS) ■ 

Once the last (finite) term is dropped from the whole amplitude, event distributions 
produced with SANC and PHOTOS become almost identical; this is shown in fig. |21 where 
comparisons of the PHOTOS (truncated SANC) with complete matrix elements are visible 
on the left-hand (right-hand) part of the figure. This is very encouraging and can be used 
as a starting point for improving PHOTOS, with the help of some kind of correction weight. 
We have indeed found that analysing the difference between the complete and truncated 
matrix elements of SANC such a correction weight 5 can be found. One can replace the 
exact matrix element of SANC with: 



E Kw^, Q,P»,Pi) ia) + MZ M (h, Q,p„, Pl )®\ 2 (1 + 6) , (7) 



where 
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E ct aa,a ; |M A ^ Ai (/c,Q, Pl/ , Pz ) (a) + M A ^ ;A! (/c,Q, M ) (fe) | 2 



-8&o(l — Pi cos i 



mfE, 



M^(l-mf/M^)(4-m'i/M^) 

, , ( ) 

+ M3 / (l-mf/M^)(4-m2/M^) ' 



and k Q denotes the photon energy, while mi and Ei denote respectively the mass and 
the energy of the fermion I. Finally M\y denotes the mass of the W boson and Pi = 
a/ 1 — mf I Ef , 6 = Z(pi,k) (all kinematical variables defined in the W rest frame). The 
differences between the results obtained from such approximated matrix elements and 
those from the complete matrix elements are rather small (see fig. 
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Figure 2: Comparisons (ratios) of the PHOTOS, truncated SANC and complete SANC pre- 
dictions for the W decay. Ratio of the [i~v acollinearity distribution from PHOTOS and 
complete SANC (left-hand side) and ratio of the \i~v acollinearity distribution from trun- 
cated SANC and complete SANC (right-hand side) are given. The dominant contribution is 
of non-leadina nature, the observable D. 



1.4r 
1.3: 
1.2 = 
1.1 : 



o 1 

re 
DC 



0.9- 
0.8^ 
0.7^ 
0.6^ 



W" -> \l v (y) 



Ratio = SANC approx /SANC exact 



-0.5 



0.5 



Figure 3: Comparison (ratio) of the predictions for the W decay from SANC truncated and 
corrected with the function 5 and from the complete SANC. Ratio of the \fv acollinearity 
distribution is given. The dominant contribution is of non-leading nature, the observable 
D. 



4 Numerical results from the improved PHOTOS 

Once our formula (jSJ) was expressed in a simple form, we could rather easily introduce it 
into PHOTOS as an additional correcting weight, to be activated in the case of leptonic W 
decays only 

As we can see from figs. HJ and |3] level of the agreement between the new version of 
PHOTOS and the complete-matrix-elements calculation of SANC is for all distributions of 
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the types A, B, C and D at the level of better than 10%. This is true, even in the case of 
the observable D, where the differences, without the correcting weights, were up to 40%. 

Also, as one can see in fig. |5J comparisons between complete and corrected-truncated 
versions of SANC give similar patterns of residual differences, as in case of PHOTOS com- 
parisons with the complete matrix elements SANC. This indicates that the approximation 
used in evaluating formula (jBJ) is at the origin of at least some of the differences between 
corrected PHOTOS and full-matrix-element calculation, rather than e.g. some technical 
problems related to its Monte Carlo algorithm. 
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Figure 4: Comparisons (ratios) of the complete SANC and corrected PHOTOS predictions for 
the W decay. Observables A and B: ratios of the photon energy (left-hand side) and muon 
energy (right-hand side) distributions from the two programs. The dominant contribution 
is of leading-log (collinear) nature. In the lower part of the plots similar comparisons for 
the complete SANC and truncated-corrected with 5 SANC predictions are given. 

With the correcting weight (J8J) included, the comparisons of PHOTOS with the "matrix- 
elements" type calculations improved significantly. The differences were not only largely 
removed, but also explained as being due to part of the W-W-j interaction explicitly 
missing in older versions of PHOTOS. This indicates that, in future, algorithms of this type 
can be improved even further, e.g. using exponentiation techniques. On the other hand, 
differences between lower and upper plots of fig. 0] point toward the limits of the method. 

At LEP2, the production and decay of W pairs is now being combined for all four 
LEP2 experiments, and uncertainties due to bremsstrahlung in W decay are important. 
PHOTOS is part of the main programs (see e.g. YFSWW3 '9\) used in the LEP2 analysis 
of the VT-pair data. The present paper provides a means of estimating the size of the 
uncertainties from PHOTOS. The correction weight can be switched on and off for that 
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Figure 5: Comparisons (ratios) of the complete SANC and corrected PHOTOS predictions for 
the W decay. Observables C and D: ratios of the photon angle with respect to fi~ (left- 
hand side) and pTv acollinearity (right-hand side) distributions from the two programs. 
The dominant contribution is of infrared non-leading-log nature for the left-hand side plot, 
and non-infrared non-leading-log nature for the right-hand side one. In the lower part of 
the plots similar comparisons for the complete SANC and truncated-corrected with 5 SANC 
predictions are given. 



purpose. We hope that, for LEP2 experiments, the updated version of PHOTOS will turn 
out to be sufficient. However, for higher precision, more sophisticated solutions, e.g. as 
those described in ref. JT0], will have to be used. Finally, let us mention that the W channel 
is of importance for some studies of the LHC Higgs-discovery potential as well pHE]. 



5 Summary 

We have successfully tested SANC versus PHOTOS in the case of the W decay. We have first 
checked that once non-leading contributions of SANC related to the W charge are dropped, 
SANC tends to agree with the results provided by the old version of PHOTOS Monte Carlo. 
These comparisons convinced us that the deficiencies of PHOTOS are indeed due to missing 
non-leading terms. The missing terms were then studied analytically within SANC and 
an approximating correcting weight, for leptonic W decays, was found. Once the weight 
was installed to the new version of PHOTOS, agreement with respect to complete results 
of SANC was found. With new correcting weight, the PHOTOS predictions are within 5% 
(instead of 7% without correcting weight) for the end parts of the spectra affected by 
the leading-log corrections, but within 5% (instead of 20%) for the angular part of the 
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distributions, where the infrared-induced logarithms dominate. The agreement is now also 
at the 10% level in the phase-space regions where only non-leading corrections contribute 
to the matrix elements. This is a significant improvement with respect to old version of 
PHOTOS, without correcting weight, where differences were up to 40%. 
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